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Abstract 

In this paper, we study the hypocoercivity for a class of linear kinetic equations 
with both transport and degenerately dissipative terms. As concrete examples, the 
relaxation operator, Fokker-Planck operator and linearized Boltzmann operator are 
considered. By constructing equivalent temporal energy functionals, time rates of 
the solution approaching equilibrium in some Hilbert spaces are obtained when the 
spatial domain takes the whole space or torus and when there is a confining force 
or not. The main tool of the proof is the macro- micro decomposition combined 
with Kawashima's argument on dissipation of the hyperbolic-parabolic system. 
Finally, a Korn-type inequality with probability measure is provided to deal with 
dissipation of momentum components. 
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1 Introduction 

1.1 Problems. In this paper, we consider the linear kinetic equation with both transport and 
dissipative terms in the form of 

dtU + Ju = lu. (1.1) 

Here, the unknown is u ^ u {t, X, with t > 0, a; e O = M'' or T'', C e R''. d>l denotes the 
space dimension. T = ^ • V^, — VxV ■ is a transport operator. V — V{x) is a given stationary 
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external forcing. L : L| — ^ L| is a linear, local, self-adjoint and non-positive operator with 
ker L ^ {0}. Moreover, L is degenerately dissipative in the sense that there is a constant Al > 
such that 

/ vLudi<-\^\\{\~-9}ufy^ (1.2) 

where Tij C L| is a Hilbert space with norm || • Ht^^, I is an identity operator and P is an 
orthogonal velocity projection operator from L| to ker L. Given initial data u(0, x, ^) = Mo(a;, 
the solution of (jl.ip formally takes the form of 

u(t) = e*^uo, B=:L-T. 

The goal of this paper is to study time-decay rates of e*^uo in some Hilbert space under some 
conditions on V and uq as time tends to infinity. 

1.2 Literature. The time rate of convergence of solutions to equilibrium is an important issue 
in the study of evolution equations. For the linear kinetic equation in the form of (jl.ip . 
the main difficulty lies in the fact that the linearized operator is degenerate in a nontrivial 
finite dimensional space since conservation laws exist for general physical models. However, 
the interaction between the transport part and the degenerate dissipative part can lead to 
convergence to equilibrium. This property is called hypocoercivity [27j . 

There have been several well-established methods to study the rate of convergence for the 
kinetic models such as the relaxation equation, Fokker-Planck equation, Boltzmann equation 
and Landau equation. Since the literature on this subject is so huge, we only mention some of 
them which are related to the study of this paper. The non-constructive method by spectral 
analysis to obtain the exponential rates for the Boltzmann equation with hard potentials on 
torus was firstly provided by Ukai [25j . The recent refinement of results of |25j can be found 
in [26] . Energy method of the Boltzmann equation was developed by Liu-Yu [21] , Liu- Yang- Yu 
[20] and Guo [131 [H] . Energy method combined with additional techniques such as velocity 
weight estimates or spectral analysis [221 [M] [S] was also applied to obtain time decay rates in 
the framework of perturbations. Another powerful tool is entropy method which works in the 
non-perturbation framework. By using this method, Desvillettes-Villani [^ obtained firstly the 
almost exponential rate of convergence of solutions to the Boltzmann equation on torus with 
soft potentials for large initial data under the additional regularity conditions. Concerning 
the Fokker-Planck equation, see [T] and references therein. In addition, on the basis of the 
spectral analysis, the hypoelliptic theory of the Fokker-Planck equation or relaxation Boltzmann 
equation was developed by Herau-Nier 15^ and Herau [14 ] I16 j . 

Recently, some general theory on hypocoercivity was provided in [22l |4l [27] . By construct- 
ing some proper Lyapunov functional defined over the Hilbert space, Mouhot-Neumann [22] 
obtained the exponential rates of convergence in i/^-norm for some kinetic models with general 
structures in the case of torus. An extension of [22] to models in the presence of a confining 
potential force was given by Dolbeault-Mouhot-Schmeiser [5 , where L^-norm was considered. 
Villani [27_ also gave a systematic study of the Hypocoercivity theory. The result on the 
Fokker-Planck in this paper is inspired by [27] . 

1.3 Notations. Define differential operators X^, (1 < i < d) by 

XiU = e 2 dxi{e 2 u), Y^-u = e ^ %(e * u). (1.3) 
Note that X.^ , Y^ are equivalent with 

X. = + Y,u = ie. + %- 

Define a norm || • by 

Mh^ = hll' + l|X«|P + ||Yu|l2 
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for u = u{x,^). Here and in the sequel, || • || means L^-norm over M.^ x M|, X = (Xi, X2, • • • , X^), 
Y=(Yi,Y2,--- ,Yd), and 

llX^f =^||X,^f, IIYuf = ^||YHP. 

i=l 1=1 

For simplicity, when a function under consideration is independent of velocity variable, || • || 
is also used to denote L^-norm over and 9, = d^t without any confusion. For the inner 
products over and iy|, we use 

{g,h)= [ f{x)g{x)dx, {g,h)= [ giOHOd^, 
respectively. For a function w = w{^), define norms | • |i„ and || • ||i„ by 

\9\l= [ NM\'+H0\9i0\'d^, 9 = 9(0, 
Il9ll^= // \Vig{x,0\^+wiO\gix,0\^dxd^, g = g{x,0. 
For g > 1 and f2 = K'*, define 

Z, = LliLl) = L\Rl;L'^{mi)), \\g\\z, = 01"^^;) d^j • 

Denote two functions M, A4 by 

M{0 = (27r)-'^/2e-l^l'/^ A^(x,0 = e-^(^'M(^), 

where M is a normalized Maxwellian. C denotes some positive (generally large) constant and A 
denotes some positive (generally small) constant, where both C and A may take different values 

in different places. A ^ B means Ai^l < B < A2A for two generic constants Ai > and A2 > 0. 
For an integrable function (7 ; f2 ^ M, its Fourier transform g is defined by 

. d 

g{k) = / e~''''''g{x)dx, x - k =:'^ Xjkj, 
Jn .^^ 

for A; e M'' if O = while kGZ'^ when n = T'^, where i = e C is the imaginary unit. For 

two complex numbers ci,C2 € C, (ci | C2) = ci • C2 denotes the dot product over the complex 
field, where C2 is the complex conjugate of C2. 

1.4 Models. To the end, we shall consider three types of degenerately dissipative operators for 
L as follows. 

Model 1. L is the linear relaxation operator 

Lu = -{I - Po}u, (1.4) 
where Pq : -L| ^ span{M^/^} is an orthogonal velocity projection operator given by 

Pou = aM^/^, a = a''=: (M^/^ u). (1.5) 

Notice that 

KerL = span{Mi/2} 
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and — L satisfies the identity 

- [ uLud^= [ |{l- PojMpdC. 

Moreover, if 

then / cquivalently satisfies the original relaxation model: 

dtf + T/ - M I fdi - f. 



Model 2. L is the linear self-adjoint Fokker-Planck operator 
It is "well-known that 

Ker L = spanjAfi/^i. 

and — L satisfies the coercivity 

- / ulud^>\{\-Po}u\l, (1.7) 

where — ly (^) = 1 + |^|2 and Po is still defined by JUS]). Instead of dUT]), it IS more conve- 
nient to use another equivalent coercivity inequality as used in [7j. In fact, define P : — > 
span{Mi/^^Mi/2} as 

Pu = Pqu® Piu, 

Pou = aAf a = a" = {M^^"^, u), 
Piw = 6-CAfi/^ 6= fe" = (^A/i/^u). 
Then, one can compute 

lu = L{\ - P}u + LPu = l{\ - P}u + LPiu = L{l - P}u - P^u. 
Thus, it holds that 

- / uLud£,> X\{\-P}u\l + \b\^. 

Similarly as before, if let 

f = M+ M^/'^u 
then / cquivalently satisfies the linear Fokker-Planck equation: 

5J + T/-V5.(V^/ + e/). 

Model 3. L is the linearized Boltzmann operator 

\-u^j^mM,M'^'u) + Q{M'/'u,M)], (1.8) 
where Q is the so-called bilinear collision operator defined by 
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for / = /(^) and g = g{^). Here, although the hard-sphere collision kernel in Q is supposed, all 
results of this paper still hold in the case of hard potentials and Maxwell molecules. For L, it 
is also well-known ^2J that dim ker L = d -I- 2, 

kerL = span{Afi/2,eiM'/',- • ' AdM'^^ ICpM^^'}, (1-9) 
and — L satisfies the coercivity 

-[ uLud£_>x[ - P}u\^d^, 

where we still used i'(^) to denote the collision frequency defined by 

HO =11 \{^-^.)-uj\Mdujd^,, (1.10) 

and also used P : L| ^ ker L to denote the orthogonal velocity projection operator. For 
convenience, corresponding to the d + 2 dimensional space (|1.9p . P is written as 

' P^i = {a + fe•e + c(|ep-d)}Ml/^ 

a — a 

' b = b' 

Therefore, a, 6, c mean mass, momentum and temperature components of macroscopic part Pu. 
If let / = + M.^/'^u then / satisfies the linear Boltzmann equation: 

at/ + T/ = Q(M,/) + Q(/,M). (1.12) 

1.5 Main results. Let us state them in two cases which will be proved in terms of different 
analytical tools. 

Theorem 1.1 (case of no force). Consider (|l.ip where Q, = or , d > 1, y = 0, and 

L is one of the linear relaxation operator, linear Fokker-Planck operator and linearized Boltz- 
mann operator as defined in Model 1, Model 2 and Model 3, respectively. Let e*^ito denote the 
corresponding solution for initial data uq = uo(x,^). 

Case (n = R'^). Let h = h{t, x, satisfy 

h{t,x,^) ±kcTL, yt>0,xeR'^. (1.13) 

Let I < q < 2. Then, for any a, a' with a' < a and m — \a — a'\, there is a constant C such 
that 

\\dSe'^uo\\<C{l + t)-^^'-i\\dS'uo\\z, + \\d^uo\\) (1.14) 

and 

\\d^ f e~^'-^^^h{s)dsr 
Jo 

< C f\l + t - s)-^^^- {}\w-^'^d^'h{s)\\\ + \\w-'/^d^h{s)f)ds (1.15) 

for any i > 0. Here, w = w(£,) is defined by 

( 1 for Model 1, 

w = u)(C) = < 1 + ICP for Model 2, 
i 1 + id for Model 3, 
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0,71(1 the index (Jam 

of the algebraic rate is defined by 

d f 1 1 \ m 



Case (Ti = T'^j. Suppose 

^iOMx,Odxd^:^0, V?AekerL. (1.16) 



Then, there are constants C and A > such that 

||e*^uo|| <Ce-^*||uo|| (1.17) 

for any t>0. 

Theorem 1.2 (case of confining force). Consider ()l.ip where ft — 'R'^, d > 1, V = V{x) is a 
confining force with 



L is one of the linear relaxation operator, linear Fokker-Planck operator and linearized Boltz- 
mann operator as defined in Model 1, Model 2 and Model 3, respectively. Let e*^uo denote the 
corresponding solution for initial data uq — uq{x,^). The following additional conditions on V 
and Uq are supposed to hold. 

Case of Model 1: = M! - |ln(27r), and 



[ [ M^/^uodxd£, = 0. 



Case of Model 2: 



^\^xV\^ -^A^V ^ (X as\x\-^(x, 
|V2F|2<<5|V,y|2 + Q, V<5>0, 

M^^^uodxd^ = 0. 

Case of Model 5; F = - f ln(27r), d>3, and 

(1, x,tx-txx^, \^\-')M'^^uodxd^ = 0, 
where {x x ^),y = Xi^j - Xj£,i, l<i,j <d. 

Then, under the above assumptions, there are constants C and A > such that 

\\e'^n^<Ce'^'\M\n^ (1.18) 

for any t > 0. 

1.6 Strategy of proof . The main idea is the macro-micro decomposition combined with Kawashima's 
argument on dissipation of the hyperboHc-paraboHc system. In fact, suppose that kerL is 
spanned by an orthogonal set 
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and the corresponding orthogonal velocity projection is denoted by 

n 

Pu = ^a^(^,a;)V'^(C)■ 
^=o 

The total energy dissipation rate corresponding to certain temporal energy functional of u = 
{I — P}u + Pu can be recovered as follows: 

Step 1. Starting from the equation (jl.ip . one can make energy estimates to obtain the dissipation 
of the microscopic part {I — P}u on the basis of the coercivity property (|1.2p of the 
linearized operator L. 

Step 2. One can derive a fluid-type one-order hyperbolic system of ai{t,x) {0 < i < n) coupled 
with { I — P}u which are actually of the moment equations in terms of the above orthogonal 
set A and some high-order moment functions. See (|2.5|) . (I2.18|) . (I4.3|l and ()4.17p for 
models under consideration. By applying Kawashima's argument on dissipation of the 
mixed hyperbolic-parabolic system 18J, one can obtain the dissipation of the macroscopic 
part Pu or equivalently ai(t, x) (0 < i < n) on the basis of the fluid- type system. 

Step 3. Combining estimates in Step 1 and Step 2, one can obtain a properly defined temporal 
Lyapunov functional which is not only equivalent with the desired total energy functional 
but also captures the total energy dissipation rate. 

Analytical tools are the Fourier transform for the case when there is no forcing and the 
direct energy estimates otherwise. When there is a potential force, we also need the Poincare 
inequality and Korn-type inequality which will be provided at the last section. 

2 Relaxation model 

In this section we prove Theorem 11.11 and Theorem ll.2l for Model 1. 

2.1 Case when n = R'^ andV ^ 0. To prove (|1.14p and p.lSp , we consider the Cauchy problem 

(2.1) 

u(0, X, = uoix, 0, xeR'^,^£ R'', 
where L is the relaxation operator given in (|1.4[) . and as in (|1.13p . h = h(t,x,^) satisfies 

h{t,x,0 -LkeiL, yt>0,xeR'^. (2.2) 
Notice that the solution to (|2.1|1 can be written as 

u(i) = e*^uo+ f e~'^*-'^^h{s)ds, (2.3) 

where B = • V^; L. 

We shall use the method of Fourier transform to deal with the time-decay of the solution u 
given by (|2.3p in a unifying manner so that (|1.14p and (jl.lSp follow from the case when h = 
and uq = 0, respectively. In fact, the Fourier transform of (|2.ip i shows 

dfii + i^ • ku = Lu + h, 

which after multiplying u and taking velocity integration and the real part, deduces 

ldt\\u\\h + \\{\-Po}u\\h^Re f {h\u)d^. 
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Due to 



(h I u)dC = / {h\{\- Po}u)dC < -\\h\\^ + -\\{\ - Po}u||i.. 



Then, it follows that 



5*11*^1112 + ||{l-Po}^^|li2 < \\h\\i., (2.4) 



which is the first estimate on the basis of the dissipative property of L. Next, we estimate 
PoM = .^yith a = a" = (Afi/^w). Recall 6 = 6" = (^Mi/^u). Then, from ^J^ji, a and 

b satisfy the fluid-type system 



(2.5) 



dta + V^-b = 0, 

dtb + v,a + V, • r({i - Po}u) = -b + /i), 

where T = (Tij)^y^ii is the moment function defined by 

T,,{g)^m,-l)M'^^g), l<t,j<d. (2.6) 
Notice by the definition of Pq that 

r({l - Po}u) = {^(g> {I - Po}u). 



Here, we used p.6p as the definition of T only for the convenience of the later proof for Model 
3. Taking further the Fourier transform of (12. 5|) gives 



( dtd + ik-b^O, 

\ dtb + ifca + ir({l - Poju) ■k = -b+ {^M^/^,h) 

where 

d 

(r({i-Po}^)-fc), = ^r,,({i-Po}u)%. 

For the above fluid-type system, from computations 

|fcp|ap = (ifca I ika) = (ifca | -dtb ~ iT{{\ - Po}u) ■k-b+ {^M^^'^,h)) 
= -dtiia I 6) + {ikdtd \ 6) 

+ (ifca I -ir({l - Po}u) ■k-b+ {^M^/\h)), 

one has 

5tRe(ia | b) + Ifcplap 
< \k ■ 6|2 + e\k\^\d\' + C,(l + |fcni|{l - Po}u\\l. + CMlp 

where e > is arbitrary. Then, taking e > small and dividing it by 1 -I- |fcp yield 

which is the second estimate based on the Kawashima's argument on the dissipation of the 
hyperbolic-parabolic system. 

Now, for i > 0, fc e R"*, define 

E{u)^\\u\\^+.Re^^ (2.8) 
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with a small constant k > to be determined later. One can let k > be small such that 
E{u) ^ 11^1112 • Taking k > further small, the linear combination of ()2.4p and (j2.7p gives 

dtE{u) + j^^Eiu)<C\\h\\^, 

which with the help of Gronwall's inequality implies 

E{u{t,k)) <e'^w^'E{uo{k)) + C f e~^^'*~'^||/i(s, fc)||2,ds. 

Hence, (|1.14p and (jl.lSp follows from the above estimate by using the standard procedure as 
in [18l [TOl [26] , and the details of the rest proof are omitted for simplicity. This completes the 
proof of prT4|) and (fLT5| for Model 1 in Theorem [TTT] 

2.2 Case when $7 = T'' and V = 0. The goal in this case is to prove (|1.17p . It actually can 
be achieved by a little modification in the case of the whole space. In fact, let's consider the 
Cauchy problem p.ip with the spatial domain M"^ replaced by T"* and h = 0. Then, ()2.4p and 
(P?7P with h = still hold for T'*. a and b still satisfy the system (^3]) . Notice that one 
has the mass conservation 

d f 

— I a{t, x)dx = 0. 
dt Jfd 

Since initially 

M^/'^uo{x,£,)dxdi ^ I aQ{x)dx = 
which corresponds to the assumption p.l6p . then 



a{t, x)dx 

fd 



= 0. 

t>o 



Thus, a(i, 0) = for any t > 0, which yields 



^^p[^|a|' > ^laP for any t > 0, fc e Z"^. 



Therefore, (12.71) is modified as 



By using the same definition of E{u) as in (|2.8p . it holds 

dtE{u) + \E{u) < 0, 

which implies ()1.17p from Gronwall's inequality and /c-integration. This completes the proof of 
pTTT]) for Model 1 in Theorem O 

2.3 Case when fl — R'* and V is confining. In particular, let V = | ln(27r). Take uq with 

[ [ M^^^uodxd^ = 0. (2.9) 

Let u{t) — e*^uo be the solution to the Cauchy problem 

(dtU + Ju = Lu, i > 0,x G R'^,f € M'', 

i (2-10) 

[u{0,x,^) = uo{x,0, xeR'^,^ewK 
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Next, we make energy estimates on u. For zero-order, it is straightforward to get 

1 d 



2 dr 



u\ 



2 



||{l-PoMr = 0. (2.11) 



For first-order, instead of directly estimating x and ^ derivatives, we use X and Y differentiations. 
Take 1 < i < d. Applying and to (|2.10p i . one has 

dtX.u -h TX,u - LX,M = [T, X,]u ~ [L, Xj]u, (2.12) 
dt^iu + TY,u - LY,w = [T, Y,]u - [L, Y,]u, (2.13) 

where [•, •] denotes commutator of two operators. Since L is local in x and X^ is only involved 
in spatial derivative and multiplier, then [L,Xi] = 0. Due to the fact that PoY^ — YiPg = 0, 

[L,Y,] = [l-Po,Y,] =-[Po,Y,] =0. 

For commutators containing T, the further computations show that 

[T,X,;]u = [e • V„X,]7/- [V,t/ • V4,X,]U 

= \i ■ VAVu + V,d,V ■ V^u =. V,d,V ■ Yu, (2.14) 

and 



= -diU - 

Thus, it follows from ((TT^ and (^1^ that 



= -diU - ^diVu = -X,u. (2.15) 



ld_ " 

2dr 



Xuf + \\{\ - Po}Xw||2 =Yj d^djVXiuYjudxd^ 



ij = l 

J^IIYuf + ||{l - Po}Yu\\' = - E // S,,X,uY,udxd^, 

where (5y is the Kronecker delta. Noticing didjV = Sij by the definition of V in the considered 
case, one has 

i|(||X«f + llY^f ) + |1{I - Po}Xu|p + |1{I - Po}Y^f = 0. (2.16) 

Since PqY ~ 0, the rest is to obtain the dissipation rate corresponding to PqXm — XPqw and 
Pqw, or equivalently Xa and a. We shall again turn to the fluid-type system satisfied by a and 
b. 

Notice that when V is nontrivial, similarly as before, from (|2.10p i . a and b satisfy 



5t(ae-^) + V, -(foe-^) =0, 

dtibe-^) + V^ae'^) + ^^Vae~^ -f • [r({l - Po}u)e-^] = -fee"* 
Equivalently, the above system can be rewritten as 

r 9ta - X* • 6 = 0, 

\ dtb + Xa - X* • r({l - Po}u) + 6 = 0, 
where X* is the adjoint operator of X given by 

V(a:) _ V(,x) 1 

XiU = -e 2 dx^{e 2 u) = {-d^,V - d^,)u, 1 < i < d. 



(2.17) 



(2.18) 
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e 2 a{0, x)dx = 0, 



= 0. (2.19) 



it follows from the mass conservation (|2.17|) i that 



e 2 a{t,x)dx 



From (|2.18p 9. one can compute 

|jXa||2 = (Xa, Xa) = (Xa, -3(6 + X* • r({l - Po}u) ~ b) 

= -^(Xa, b) + [Xdta, b) + (Xa, X* • r({l - Po}u) - b) 

where it further holds from (|2.18|) i that 

{Xdta, b) = {dta, X* ■b) = (X* ■ b,X* ■ b) = \\X* ■ bf. 

Then, it follows that 

^{Xa, b) + A||Xa||2 < ||X* • bf + C(||X*{I - Po}u||2 + \\bf). (2.20) 
at 

Lemma 2.1. As long as there is a constant C such that 

lA^Fp < C(|V,Fp + 1) (2.21) 
for all X G K*^, there is some constant C such that 

\X*g\^dx <C I (IX5P + \g\^)dx (2.22) 

for g = g{x). 

Proof. From integration by parts, 

j |X*5pdx= / \Xg\'^dx+ [ \g\^A^Vdx. (2.23) 

jR'i JRii JR'i 

Under the assumption (|2.21[) . 

\g\^A^Vdx< [ \gme\A,V\^ + ^)dx 
.had 4e 



< / |g|2(C7e|V,y|2 + Ce+l)dx 

= 4Ce/ |<?p(i|V,t/p-iA,T/)dx + 2Ce / \g\''A,Vdx 



+ (Ce+l) / lal'dx 



for arbitrary e > 0. Taking e > small, one has 

[ \g\^A.^Vdx<C ( \Xg\^dx + C [ \g\^dx. 
Then, (|2.22|) follows by plugging the above inequality into (|2.23p . □ 
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So, by applying Lemma [2TT] to (|2.20p . one has 



-{Xa,b) + A||Xaf < C(||{l - Poj^uf + \\{\ - Po}"ll'). 



Furthermore, due to (|2.19p . one has Poincare inequahty 



llXaf >A||a|p 



from Proposition 1 5. II Then, it foUows that 



^(Xa, b) + XiWXaf + \\af) < C(|l{l - Po}Xu\\' + \\{\ - Po}u\\'). 



(2.24) 



Now, let us define a temporal functional 



£{u{t)) = + ||Xu||2 + ||Yu|p + K(Xa, b) 



with a small constant k > to be determined. Firstly, k > is chosen small such that 



K > is further small enough such that the linear combination of ()2.1ip , (|2.16p and ()2.24p gives 



By Gronwall's inequality, this proves (|1.18p and hence completes the proof of Theorem 11.21 for 
Model 1. 

3 Fokker- Planck equation 

In this section we prove Theorem 11.11 and Theorem 11.21 for Model 2. The proof in the case of 
no forcing can be carried out in the completely same way as for Model 1. When a stationary 
potential forcing is present, the proof for Model 1 can be refined to yield the exponential 
time-decay rate in even for a class of general potential functions essentially because the 
Fokker-Planck operator enjoys the velocity regularity. This actually has been studied in detail 
by Villani [27], and the hypoelliptic theory is founded by Herau [Ml [16]. Here, we shall give 
another proof which is based on the macro-micro decomposition and Kawashima's dissipation 
argument on the hyperbolic-parabolic system. 

3.1 Case when = M.'^ or T"^ and V — 0. We consider the Cauchy problem 

< (3-1) 

where L is the self-adjoint Fokker-Planck operator given in (|1.6p . and h = h{t,x,£^) satisfies 

a;, _L ker L, yt>0,xen. 
As before, the solution to p.ip can be written as 



where B = — ^ • + L. In this case. Theorem 1 1.1 1 and Theorem 11.21 for Model 2 can be proved 
in the same way as for Model 1. Thus, all details of the proof in this case are omitted for 
simplicity. 



S{u{t))^\\u{t)rm. 



^£{u{t)) + \£{u{t)) < 0. 
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3.2 Case when H, — M.'^ and V is confining. Suppose 

-^(")da; = 1, (3.2) 



e 

j\\7,V\^ ~^A,V ^ ^ as |a;|^oo, (3.3) 

\'^lV\<S\\/^V\^ + Cs, V(5>0. (3.4) 

Take uq with 



Let u{t) = e*^uo be the sohitio to the Cauchy problem 

j dtu + Tu = Lu, t>0,xeR'^,^eM.'^, 



(3.5) 



(3.6) 



Next, we follow the same line of proof as for Model 1. It is straightforward to get from 
that 

~\\ur + X\\{\-P,}u\\l<0, (3.7) 

where — i^(^) = 1 + Recall the definition (|1.3p of the operator Y. Observe that since 
Lu = i(2d- ICpjw, then [L, Y] = Y. Thus, similar to obtain (|2rT2)) - ((2?T3)) with identities 

([214)1 - (|2T^ . one has 

dtX,u + TX,u - LX,it = V^d^V ■ Yu, 
dtYiU + TYiU — LYiU — ~XiU — YiU, 



which imply 



1 d 
2di 



\\Xuf + \\\{\ - Po}Xu||2 < V // d^djVX^uYjudxdC, 



(3.8) 



].^\\Yuf + \\\{\-P^}Yu\\l + \\Yuf<~y^[[ X,uY,udxdt (3.9) 



The r.h.s. of (|3.8p is bounded as 

d 



J2 d^djVX.uYjudxd^ 

= V // d,djV{\~Po}X,uY,{\-Po}udxd( 

= V // d^djVY*{\ - Po}X,;u{l - Pojudxd^ 

< e||{l - Po}X7.||^ + - // \^lvf\{\ - Po}ufdxd^ 

for an arbitrary constant e > to be chosen later, where we used YPgu = and Y* is the 
adjoint operator of Y given by 

Y*u = -e—d^^{e~ — u) = (-^j - 
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Furthermore, due to the assumption (|3.4p on V, 
II \^lV\^\{\-9o}u\^dxd^ 

' // \V.V\^\{\~P^}u\^dxdi + Cs 1 1 \{\~9o}u\^dxd£, 



< S 



4<5 / / (^|V,Fp - ^A,V)\{\ - P„}u\^dxd^ 



+25 1 1 A^V\{\-Po}u\^dxd£, + Cs 1 1 \{\ - Po}u\^dxdi, 
which by smallness of 5 > impHes 

\VlV\^\{\ ~ Po}u\^dxd^ < CS\\{\ - Po}Xu\\^ + Cs\\{\ - Pojuf. 



So, from the above estimates, (|3.8p is bounded by 

^l\\Xuf + X\\{\-Po}Xu\\l 

< + ^) ll{i - PolX^^II^ + ^||{l - PoM'. 

Since e > and 5 > can be arbitrarily smaU, it foUows that 
1 d 



2 ^JXuf + X\\{\ - Po}Xu\\l < C\\{\ ~ PoMp. (3.10) 
Noticing 

(X,u,Y,u) = ({I-Po}X,u,Y,m), 

one has from (|3.9p that 

^|l|Y«f + A||{l - Po}yu\\l < C\\{\ - Po}Xur. (3.11) 

FinaUy, we turn to the estimate on PoXu and Pqu. Similar to obtain (|2.18p . (|3.6P i also gives 
the same fluid-type system: 



( dta^X* ■b = 0, 

1 dtb + Xa-X* ■ r({l - Po}u) + 6 = 0, 



with 

0. 



e 2 a(t, x)dx 



t>o 



again due to the mass conservation and initial condition p.Sp . Applying the same argument 
as before and then using Poincare inequality in Proposition 15 . 1 1 bv assumptions (|3.2p - (|3.3p . one 
has 

^(Xa,6) + A(||Xa||2 + ||a||2) < C(||{l - Po}Xuf + \\{\ - Po}uf). (3.12) 
at 

Now, let us define a temporal functional 

£{u{t)) = \\uf + Ki\\Xu\f + K2\\yuf + K3iXa,b) 
with constants > (1 < i < 3). One can choose 

< K3 < K2 < < 1 
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such that 

and further the linear combination of (|3.7p . (|3.10p . (|3.1ip and (|3.12p gives 

^^£{u{t)) + X£{u{t)) < 0. 

By Gronwall's inequality, this proves (|1.18|) and hence completes the proof of Theorem 11.21 for 
Model 2. 

4 Boltzmann equation 

In this section we prove Theorem 11.11 and Theorem 11.21 for Model 3. Although this can be 
done along the same line as for Model 1 and Model 2, it is a little more complicated for 
Model 3 since ker L = d + 2 so that the linearized Boltzmann operator L is degenerate over a 
space with higher dimensions. The key idea is still based on the macro-micro decomposition 
and Kawashima's dissipation argument on the hyperbolic-parabolic system. The additional 
difhculty in the presence of confining forces lies in verifying the Korn-type inequality to obtain 
the dissipation of momentum component in the fluid part, which is left to the next section. 

4.1 Case when il — R'^ and V — 0. We consider the Cauchy problem 

dtu + £, -Wxu ^ Lu + h, t > 0.x eR"^.^ eM.'^, 

(4.1) 

u(0, X, = uaix, 0, X e R'*, C e M'', 

where L is the linearized Boltzmann collision operator given in p.Sp . and as in (|1.13p . h — 
h{t, X, £,) satisfies 

h{t,x,^) J-keiL, yt>0,xeR'^. 
As before, the solution to (|4.ip can be written as 

u(i) = e*^uo+ / e-^'~'^^h{s)ds, 
Jq 

where B = — ^ • Vx + L- Both parts in the above u{t) will be estimated in a unifying way by 
using the method of Fourier transform similarly before. 
Firstly, from (|4.ip i . it is straightforward to obtain 

U\\u\\h + A // - P}u\'dxdi < C\W-'/'h\\l. (4.2) 

where ~ 1 + defined by (|1.10p is the collision frequency for the case of hard sphere 
model [To], and the orthogonal velocity projection operator P : L| ^ ker L is described by 
()1.11|) . Next, we devote ourselves to the estimate on Pu or equivalently (a,5, c). In fact, by 
taking ^d{d -I- 5) + 2 number of velocity moments 

Afi/2,e»Mi/2, (lep - d)M'/\ (e,o - i)M'/', (lep - d - 2)e.Mi/2 

with 1 <i,j < d ior the equation (|4.ip i . one has the fluid- type system 
' dta + Wx-b = 0, 
dtb + V,(a + 2c) + V, • r({l - P}u) = 0, 

< atc-K^V, •5+^V,-A({l-P}u) = 0, (4.3) 

dt[r,,{{\ - P}u) + 2c%] -I- + djb, = r„(r + h), 
^dtA,{{\-P}u) + d,c = A,{r + h), 
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where a, b, c are defined in (|l.lip . the matrix- valued function F is defined by (|2.6p . the moment 
function A ~ {■^i)i<i<d is defined by 



A.(5) = ((ICp-rf-2K.MV2,5), 



and r is denoted by 



r = V:,{l - P}m + L{l - P}u. 
The detailed derivation of (|4.3p was given in [6] and thus is omitted for simplicity. 

Lemma 4.1. There exists a functional Ei-^tiu) given by 

d 



1 



^ {ikSj + ikfh I r,j({l - P}u) + 2c4 



'«2 



d 



■ ^(ifciS I 5i) 



dtc +-,ik-b+ —\k ■ A({l - P}u) = 0, 
a 2d 

dt[T,ji{\ ~ P}u) + 2cSij] + ikibj + ikjk r,y (f + h), 
{ dtA,{{\ - Pjii) + ihc ^A,{f + h), 
for 1 < i, j < d, where 

f = -i^ ■ k{\ ~ P}u + L{\ ~ P}u. 
Step 1. Estimate c: Starting from the highest-order moment equation M.?!) -;. one has 

d d 



|A:|2|c|2 = Y^iihc I ihc) = ^(ifc.c I -atA,({l - P}^) + A,(f -f /i)) 

i=l i=l 

d d 

= -dtY^iihc I A,({l - P}u)) + Y.^ihdtc I A,({l - P}u)) 

i=l 1=1 

d 

+ Y^{lhc I ^^{f + h)). 



(4.4) 



(4.5) 



' ' i— 1 

/or constants < ^2 ^ /^i ^ 1 suc/i ^/lai t/iere are constants A > 0, C suc/i ^/?-ai 

dtReE-Uu) + + I^P + < - P}"lli| + \W-"^h\\li) (4.6) 

/loWs for any t > and fc € M''. 
Proof. Write Ij4.3p in terms of Fourier transform as 
dtCL + i/s • 6 = 0, 

dtb + i(a + 2c) + iF({l - P}u) ■ k = 0, 



(4.7) 
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Using (|4.7p q to replace dtc gives 

d 

Y,{ikAc I A,({l - P}u)) 



i=l 

d 



= X^(^ifc -b+^ik- A({l - P}u) I ifc,A,({l - P}u)) 



2d 

i—1 



<«ilt-il' + -ll{i-p)i|li,, 



2| 

where < ei < 1 is arbitrary to be chosen later. Moreover, it holds that 



^(ifc,c I A,(f + h)) < lifcHcp + c(i + i^niKi - p}u\\^ + c\w-'/'h\\l.^ 

i=l 

Then, it follows 

d 

dtRe ^(ifc,c I A,({l - P}u)) + Alfc^cl^ 

i=l 

< eilfc • 6p + -(1 + \km\ - P}u\\l. + C\W-'/^h\\l.. 

Step 2. Estimate h: Observe the identity 

d 

^ \\\hhj + i/cj6,|p = 2|/cn6|2 + 2||fc • 6f . 

On the other hand, compute from (|4.7|) zi that 

d 

^ |lifc,6j +i%&,f 



= {ikib, + ikjk I -^t[T^J{{\ - P}u) + 2c%] + (f + /i)) 

d 

= -5t ^ iikibj + ikjk I Ey ({I - P}u) + 2c%) 

d 

+ J2 iShdthj + ikjdtk I ({I - P}u) + 2c%) 



- ^ (ifcifoj + i%&j I Ey (f + /i)). 



Using (|4.7p 9 to replace 9(6 gives 

^ (ifc,5t&j + ifc,9t6, I E„({l - P}u) + 2c%) 



ii=l 

d 

= 2 ^ (9^6, I i%E,, ({l - P}u) + 2ifc,c(5,,) 

ij = l 

d d 

= 2 ^ (ifc,,(a + 2c) + i^ fc£E£,({l - P}u) \ \kjT^j{{\ - P}{t) + 2i%c4 

ij = l 1=1 

< 62|fcna|2 + ^|A:nc|2 + -|fcni{l - Pj^lli., 

£2 €2 « 
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where we used the symmetry of and < 62 < 1 is arbitrary to be chosen later. Moreover, 
it is straightforward to obtain 



(ifci^j + ikjbi I Tij{f + h)) 



ij = l 



d 

< ^ E W'^^^i + '^M? + c{i + ifcnii{i - p}u\\i. + c\\i,-^'^h\\i^ 



2 

So, coUecting the above estimates, one has 

d 



dtRe E i^kSj + ikfb, I ry({l - P}{t) + 2c(5y) + |fc|2|6|2 

y = l 

< e^lfeplap + -lA^Plcp + -|fcn|{l - P}u\\l, + C||^-V2^||i.. (4.9) 

£2 £2 ^ « 



^tep 5. Estimate a: Similarly, ([4.7^ 9 implies 

I^HaP = E(i^'^ I = E(^^'" I ~ E i^^y ({I - P}u)) 

i=l i=l 
d d 

+ ^{iha I -2ifcjC ~^ikjTij{{\ - P}u) 
1=1 j=i 

which after using (|4.7P i to replace dta on the r.h.s. gives 

d 

dtRe ECi'^^a I b^) + X\k\^\a\^ < \k ■ + C\k\^\c\^ + C\k\^\\{\ - PjuH^s. (4.10) 

Finally, define Eint{u) by (|4.5p . Then, for properly chosen constants < K2 ^ ki <C 1 
and small constants ei > 0, e2 > 0, (|4.6p follows from the linear combination of the above three 
estimates (|4.8p . (j4.9p and (I4.10p and further dividing it by 1 + This completes the proof 
of Lemma 14.11 □ 

Remark 4.1. As in ()2.7|) . p.24p and p.l2p . Lemma \4-.1\ shows that the dissipation of degener- 
ate part which is the kernel of the linearized Boltzmann equation with dimensions equal to d + 2 
can be recovered from the fluid-type moment system (|4.3p . This property was firstly observed by 
Guo \1SX [TBI, '^'^'^ ^'^■s Iciter improved in JMlBtiSj and !1L , J7 / for different purposes. Precisely, 
US^ is mainly based on an elliptic-type equation ofb derived from HA.S^ i- ()4.3P a. The aim of |^ 
is to remove time derivatives by constructing some functional which is similar to (|4.5p but takes 
more complicated form. /6/ introduced moment functions F and A to refine the form of (|4.5p . 

exactly used the same method to deal with the Vlasov-Poisson-Boltzmann system. Here, it 
should be emphasized that due to new estimates on the dissipation of b, the current method of 
proof is more general in the sense that it can be directly modified to apply to the case with a 
potential forcing, see Lemma Notice that three terms in ()4.5p are inner products of i-th 

order moment and [i + l)-th order moment for i — Q, 1,2, respectively, and also it is in the 
same spirit of Kawashima 's construction of compensation functions \10^ . 

Now, for t > 0, fc e R'^, define 

E{u) = \\u\\l2 K3ReEint{u) 



Hypocoercivity of Linear Degenerately Dissipative Kinetic Equations 



19 



where £'int(u) is defined by (|4.5p and K3 > is to be chosen. Similarly before, one can let 
K3 > be small such that E{u) ^ ||u||^2 and the linear combination of (14. 2p and (j4.6p gives 



dtE{u) 



,E{u)<C\W-'^'h\\l,^ 



l+|fc|2 

Hence, (I1.14[) and (|1.15p follows from the above estimate in the same way as before. This 
completes the proof of (|1.14p and (jl.lSp for Model 3 in Theorem 11.11 

4-2 Case when ~ T"^ and V = Q. In this case, (11.171) follows from the same argument as for 
Model 1 and Model 2. In fact, it suffices to check 



(a, 6, c)dx 



= 



t>o 



which results from the fact that it initially holds at t = by (|1.16p and the conservation laws 

d f 

— / (a, b, c)dx — 
at Jfd 

holds by ^^i-^^-^. This completes the proof of pTT7)l in Theorem O for Model 3. 

4-3 Case when f2 = K'* and V is confining. In particular, let V — | ln(27r) and d > 3. 

Take uq with 

(1, X, C, a; • x X \^f)M^/'^uodxd^ = 0, (4.11) 



where {x x f)y = XiS,j — XjS,i for 1 < i, j < d. Let u{t) — e*^uo be the solution to the Cauchy 
problem 

{dtu + lu = i.u, i > 0,x e R'*,£ G M'', 

{ (4-12) 

[ u(0,a;,^) = uo(a;,0, x CM!^ . 

Firstly, we verify that the property (jl.lip can be preserved at all time t > for the above 
Cauchy problem (|4.12p . 

Lemma 4.2. Under the assumption (|4.1ip . 

(1, X, a; • a; X |a;|^ \£,\^)M^/'^u{t, x, i)dxd£, = (4.13) 
holds for any t>Q. Particularly, for any t>Q, 



e ^a(t,x)dx 



e 2 li^(t^ x)dx 



e 2 c{t, x)dx = 0, 1 <i <d, 



e 2 a; X 6(t, x)dx — 0. 



Proof. Let ((•, •)) denote the usual inner product in ^. For simplicity, define also the weighted 
inner product ((•, •))a< by 

((u,w))a< = / / uvM^^^dxd^. 

Notice that f = A4 + M.^^'^u satisfies (I1.12p . From integration by parts, it is easy to obtain the 
following ODE system 

|((l,/)) = 0, |((a;xe,/))=0, 

= ((?,/)), j^UJ)) = -{{xJ)). 



-\x\\f% 



dt 

|((|xp, /)) = 2ix ■ /)), |((|ep, /)) = -2ix ■ t /)). 
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Equivalently, in terms of u, it holds that 



^((1,w))a4 = 0, ^{{x X ^,u))m =0, 



(4.14) 



dt 



([x ■ ^, u))m 



\x\'^,y'))M, 



j^(i\x\',u))M = 2{{x-tu))M, j^im^u))M = -2(ix-C,u) 



Ml 



where we used the fact 



and 



Define temporal functions 

yi{t) ^ {{S.i,u))M (l<«<rf), z{t) ^ {{x ■ (,,uf}M- 

Then, from (|4.14p and (|4.13p . yi{t) and z{t) satisfy the initial value problems of the linear 
second-order ODE as follows 

2/-'W+2/»W = 0, t>Q 

2/i(0) = -((a;z,Mo))x = 0, 

' +4z(t) = 0, t>0 
2(0) = Clio ))ai=0, 

{z\Q)^~i\e-\x\\u,))M=Q. 
Hence, both solutions are trivial, i.e. 

y,[t) = (1 < i < d) and z(t) = 

for any t>Q. Putting this back to the system (I4.14p implies that all inner products in (|4.14p 
vanish due to (|4.13p . This completes the proof of Lemma [4.21 □ 

Next, we proceed along the same line of proof as for Model 1. The zero-order energy 
integration of (|4.12p i gives 



1 d 
2di' 



|wf + A // J/(OI{l - P}u\'^dxd^ < 0. 
Similar to obtain p.l6p . for X and Y differentiations, one has 

J^(l|Xu|p + IIY^f ) + A // - P}Xu\'dxd^ 



(4.15) 



-l-A 



[[ i^{C)\{\-P}Yufdxd(<~y2 [[ [L,Y i\uYiudxdC 



It is easy to check Pu E ker L which implies 

[L,Y,]u = LY,u -Y,lu= IY,{\ - P}u - Y,L{I - P}u = [L, Y,]{l - P}u. 
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Since [L, Y^] is a bounded operator from i| to by [9], 



E 



[L,Y,]uY^udxd^ < e\\{a,b,c)\\^ + -\\{\ - P} 



for an arbitrary constant < e < 1 to be chosen later. Then, it follows 



+A 



- PjYul^dxd^ < e\\{a,b,c)f + -\\{\ - P}u 



(4.16) 



The rest is to deal with the dissipation rate corresponding to PXu and Pu or equivalently 
X(a, 6, c) and (a,b,c). Poincare and Korn-type inequalities given in the next section play a 
key role in this step. In fact, similar to get (|4.3p . from the equation (|4.12p i in the presence of 
potential forcing, one can obtain the following fluid-type system: 

( dta~X* - b^O, 
dtb + X(a + 2c) - X* • r({l - P}u) = 0, 

dtc+^X-b-^X* ■A{{\-P}u)^0, (4.17) 

dt[T,,{{\ - P}u) + 2c%] + X,6, + X,b, = (r), 

[dtA^{{\~P}u) + X,c^A,{r), 

where moment functions F and A are defined as in (|2.6|) and (|4.4p respectively, and for simplicity, 
we still used r to denote 

r = -^- yx{\ - P}iL + "^xV ■ VJI -P}u + L{\- P}u. 

Lemma 4.3. There exists a temporal Junctional £int{u{t)) given by 



fintKi)) =^(X,C,A,({I-P}li)) 

i=l 

d 

+Ki J2 O^^bj + Xjb,, ({I - P}u) + 2c%) 



ij=l 
d 



+K2 ^(X^a, bi 



(4.18) 



for constants < K2 ^ ki ^ 1 such that there are constants A > 0, C such that 

^AnMt)) + XV{a, b, c) < C{\\{\ - P}Xuf + \\{\- P}uf) (4.19) 
dt 

holds for any t > 0, where the dissipation rate 'D{a, 5, c) denotes 

Via,b, c) = (llXaf + \\a\f) + J^Wl' + M') + iWM' + M')- 

Proof. This follows from the same procedure as for the proof of Lemma 14.11 which is on the 
basis of Fourier transform when there is no external forcing. For completeness, we shall provide 
all details of proof. Firstly, it is straightforward to check 



\\T{r)f + \\A{r)f<C 



|V,{I - P}uf + \{\ ~ P}uWl + \VxVf)dxd^ 



< C{\\{\ - P}Xi 



\{\-p}ur). 



(4.20) 
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Step 1. Estimate on c: (|4.17P k gives 

d 

||Xc||2 = ^(X,c, -dtA,{{\ - P}u) + Ur)) 

1 d d 

= -- 5](X,c, A,({l - P}^.)) + Y.{XAc, A,({l - P}u)) 

i=l i=l 

d 

+ ^(X,c,A,(r)). (4.21) 

1=1 

Using (|4.17p ^ to replace dtc, one has 

d 

Y,{XAc,K,{{\-P}u)) 

= (-^X . 6 + ^X* . A({l - P}^), X* . A({l - P}u)) 

<ei||X.6f + -||X*{l-P}u|p 
ei 

for an arbitrary constant < ei < 1 to be chosen later. In addition, it holds 

d 

E(X.c,A,(r))<-||Xcf + qiA(r)|p. 

i=l 

By putting the above two estimates into (|4.2ip and using (|4.20p and Lemma I^TTl one has 
, d 

5^(X,c,A,({l-P}u)) + A(||Xc||2 + ||cf) 



dt 

i—1 



,2 . C ^ 



< eillX . + -(||{l - PjX^II^ + ||{l - P}uri (4.22) 
ei 



where due to Lemma l42l and Poincare inequality in Proposition l5.ll |lc|p was included in the 
dissipation rate. 

Step 2. Estimate on b: (|4.17P zt gives 

d d 

= l ij=l 

, d 

d 

+ E (X^atfe, + X,dth,T,,{{\ - P}u) + 2cdy ) 

d 

+ E(XA+X,5.,r,,(r)). (4.23) 



Hypocoercivity of Linear Degenerately Dissipative Kinetic Equations 23 
From (|4.17[) 9. dtb is replaced to obtain 

d 

J2 i^^^tbJ + X,dth,r,,{{\ - P}u) + 2c(5y) 

d 

= 2j2idA,X*r,,{{\ - P}u) + 2X*cS,,) 

ij=l 

d d 

= 2 ^(-X,(a + 2c) + ^X|r,,({l - P}w),X*ry({l - P}u) + 2X*c%) 

y"=l t=l 

< e2||Xa|p + -(||X*c||2 + ||Xc||2) + -||X*{I - P}u\\^ 

£2 £2 

for an arbitrary constant < €2 < 1 to be chosen later. The final term on the r.h.s. of (j4.23|) 
is bounded by 

d 1 ^ 

^(x,6, + x,&„r,,(r)) < - ^ \\x,b, + xM' + c\\r{r)r. 

ij=l ij = l 

Similar to Step 1, it follows from (|4.23p that 

J d d 

- ^(X,&, + X,6„r,j({l - P}u) + 2c%) + A;^(||X6,||2 + ||5,^||2) 

< e^llXaf + -dlXclp + \\cf) + -(||{l - P}X«f + ||{l - P}uf), (4.24) 

£2 £2 

where due to Lemma [4.21 Korn inequalities (|5.5p and (|5.14p were used. 

Step 3. Estimate on a: (|4.17p 9 implies 

d d 

\\Xaf = 51 (^^«' - 2X»c + - P}u)) 

i=i j=i 

id d 

= --^(X,a,6,) + ^(Xi9ta,60 

i—l i=l 
d d 

+ ^(X,a, -2X,c + ^ X*r,,({l - P}u)), 

i=i j=i 

where it further holds that 

d d 

Y,i.^^^ta,b,) = {dta,X* ■ b) = \\X* ■ b\\^ < C^dlXfeJ^ + ||;,^||2) 

i=l 1=1 

by dUTDi, and 

d d 

^(X,;a,-2X,c + ^X*r,,({l - P}u)) 

t=i j=i 

< ^\\Xaf + C\\Xcf + C\\X*{\ - P}uf. 

Then, it follows 

-|^(X,a,M + A(||Xa||2 + ||a|r) 



dt 

1=1 

d 

< Cj2m,r + \M') + CWXcf + C{\\{\ - P}Xu\\' + ||{l - P}ur), (4.25) 



i=l 
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where again due to Lemma |4.2[ Poincare inequality in Proposition 15. II was used for a. 

Finally, define £int{u{t)) as in (|4.18p . Then, for properly chosen constants < K2 <C 
Ki <C; 1 and small constants ei > 0, 62 > 0, (|4.19p follows from the linear combination of three 
inequalities (|4.22|) . (|4.24p and (I4.25p . This completes the proof of Lemma l43l □ 

Now, for t > 0, define 

£{u{t)) = \\uf + «3(llX?/f + llYuf) + K4fi„t(u(i)) 

where £i^t{u{t)) is defined by (|4.18p and constants < K4 ^ K3 <C 1 are to be chosen later. By 
letting < K4 ^ K3 ^ 1 be properly chosen, one has 



£{u[t)) ^ \\u{t)\\ 



On the other hand, by further properly choosing < K4 ^ K3 ^ 1, the linear combination of 
(I4T511 . (j4T6)l and (l4T9l) yields 

j^8{u{t)) + X£{u{t))<0, 

where e > in (|4.16p was also chosen small enough. Hence, p.lSp follows from Gronwall 
inequality. This completes the proof of Theorem 11.21 for Model 3. 

5 Appendix: Poincare and Korn-type inequalities 

Let V = V(x) be a smooth confining function over M'' with 

e-^(^)dx = L (5.1) 
We need assumptions on V: 

^IVV^P - ^AF -> 00 as 00 (5.2) 

or 

1 1 

4I VV^P - -AV - \9idjV\ ^ 00 as |a;| ^ 00. (5.3) 
ij=i 

Define operators (1 < z < d) associated with V as in (|1.3p . Throughout this section, 
differential operators V,A,V^ denotes W x, ^x,'^x foi' simplicity since all functions are those 
of spatial variable. The following Poincare inequality in Mf^ has been used in this paper. 



Proposition 5.1. Let d > I. Suppose (|5.ip and (|5.2p hold. There is a constant A > such 
that 

J2 f \X,a\^dx>xf \a\^dx (5.4) 



holds for any a : M'* ^ M G L^(K'*) satisfying 



1=1 

2(Tad 



V(x) 

e 2 adx = 0. 



For the proof of the above proposition, refer to [27] on the basis of the constructive method, 
where A can be explicitly computed. 

The following Korn-type inequality in R"* was also essentially used in the previous proof to 
obtain the dissipation of the macroscopic momentum. 
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Theorem 5.1. Let d > 2. Suppose (|5.ip and (|5.3p hold. There is a constant A > such that 

d 



V/ \X,bj+Xjb,\^dx>Xy2[ \b^\^dx 

holds for any b = 62, • • • , 6d) : M'^ ^ M'^ e (^^(R'i))^ satisfying 

f ^(^) 

/ e ~b{x)dx — 

and 



(5.5) 



e 2 — Xjbi)dx — 0, I < i =/= j < d. 

Proof. We use the contradiction argument. 

Step 1. Otherwise, for any n > 1, there is a non-zero element 6" — (6", ' ' ' ' ^d) ^ (^^(K''))'* 
with 

/ e — ^6"(x)dx = 

and 



_ V(a;) 

e ^ XX b"{x)dx = 0, 



such that 



ij = l 

W.l.g., one can suppose 



Thus, 



J2 [ \X,b-+X,b-\^dx<-J2 [ \b7\'dx. 
V / \X,b]+X,b':fdx<-. 



ij = l ' 

S'tep ^. We claim that there is a constant C independent of n > 1 and a continuous function 
w{x) only depending on V with 

inf w{x) > 1 and w{x) ^ 00 as —^ 00 

such that 



sup / |V6"|2 + \b''\^w{x)dx < C. 

n>l jRd 



In fact, one can compute 

d 



V / \X,b]+X,b^\^dx^2y2 I |X,67pdx + 2 V(X,67,X,6r) 

^■.h+h- (5.6) 

For /i, it holds that 

d r. 

h^2y] / \x,b]\''dx 
= 2^ / m\' + \b-f{\mf-ld.Av)dx 

tj — ± 

= 2 f iVfe'f + |6'f (ilVFp - iAy)dx. 
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For I2, from integration by parts, it holds that 

d d 

'2"" ' 2 



ij=l 
d 



dd.Vb], ^d.Vb-) + [d.bl d,b-) + dd.Vbl d,b-) + {d,bl Id.Vb^) 



ij = l 
d r 

lj=l 

d 



(.ld,Vb7,ld,Vb]) + id.b-,d,b]) 



-{ld,d,Vb],b:) - dd.Vd.blb'l) - d^^^,Vb-,b-) - {ld,Vd.b-,b]) 



which further can be written as 
I2 = 2 



(^VV^ • 6", ivy • 6") + (V • 5", V • &") - (Vy • 6", V • 6") 



-2(6", V^yfe") 

2|| ivy • 6" - V • - 2(6", V2y6"). 



Thus, it follows 



/i + /2 = 2 / |V6"|2 + |6"|2(i|vyp - iAy)da 
+ 11 ivy • 6" - V • 6"|p - 2(6", V2y6") 



Then, this implies 



Ii + /2 > 2 / |V6"p + rfi\\VVf - iAy - E |ai5,y|)da 



+ ||_Vy -6"- V-6"||2. 



Now, define 



wix) = max{i|Vyp - iAy - E \9^djV\, 1}. 



ij=l 



Due to (|5.3p . w{x) ^ cx) as |a;| ^ cx). Choosing R> large enough, one has 



;(x) = i|vyp-iAy-El5,:5,y| 



ij = l 



for Ixl > ii. Then 



/ r\'{j\VV\'~lAV-J2mV\)dx^ [ +[ 



<R 



|6"|^u'(x)dx- 



|6"|^w(x)dx + 



2;|<_R 



x|<_R 



> / |6"|^w;(x)dx-Cfl,y, 
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where 



Cb,v = sup w{x) - inf (i|VF|2 - l-^v- V \d,d,V\). 

\x\<R l^l<« 4 2 

Then, it follows that 



That is, 



< -Z \- Cr v < n + Cr,V- 

2n 2 



Hence, the claim follows. 



Step 3. From the Rellich-Kondrachov compactness theorem, up to a subsequence, (&")„>i 
strongly converges to some function b e (L^(M'^))'*, i.e. 



Thus, h satisfies 



and 



- fell as n oo. 

||6|| = 1, (5.8) 

X,h,j + Xjhi = 0, l<i,i <d, (5.9) 

e"^6(a;)dx = 0, (5.10) 



e 2 (xifoj — Xjbi)dx — 0, I < i ^ j < d. (5-11) 



We claim that (|5.9p . (|5.10p and (|5.1ip imply 6 = 0, which is a contradiction to 
In fact, set 

bt{x) ^ e~^bi{x), i = l,2,---,d. 

From Xibi ~ 0, 

ftfe,, = 0. 



Then, bi{x) is independent of Xi. Thus, set 6; — bi{xi), where Xi is a variable (xi, a;2, • • • , a;^) 
excluding Xi. From (|5.9p for i ^ j , one has 

a,5,(ij) + c)j6,(i0 = 0, z^j. (5.12) 

This implies that bi{xi) is linear in each (j 7^ i). Thus, 

{q: — 0,aj —0 or l.j^i} 

It suffices to prove all coefficients (7^ vanish. In what follows, fix z G {1, ■ • ■ , c?}. 
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Case 1. \a\ = 0. Due to ((5l0l) . 

Case 2. 2 < |a| < d - 1. Let I = \a\ and take 

{«i,«2, • • • C {1,2, • • • 
Consider the coefficient of the monomial 

It only appears in hj with j £ {1, 2, • • • , 12, • ■ ■ , it}- Take 

e {1,2,--- ,d}\{ii,i2, ■ • • 
Consider the following three functions 

" ■ ■ ^" C'2'^zi*^i2 ' ' ' '^ii — i'^ii-i-i ' * * ; 
^2^—1 ■ ■ ■ ^" C^Xi^X I2 ' ' ' Xi^_2Xif Xi^^-^ ~t~ • • • , 

where Ci, C2, C3 are the corresponding constant coefficients. Due to (|5.12p . 

Cl + C2 = 0, C2 + C3 = 0, C3 + Cl = 0. 
Thus, Cl = C2 = C3 = 0. The coefficient of ^'j^iXi- is zero. This proves 

C = o 

when 2 < |a| < d - 1. 

Case 3. \a\ = 1. Take j ^ i- Consider the function 

hj = h CiXi H , 

where Ci is a constant. Due to (|5.12p . hi is in the form of 

bi = ■ ■ ■ — CiXj + • • ■ . 

By using (|5.1ip . one has 

'" I „ r \ . I „-V/ 2 , 2 



= / e {xihj — Xjbi)dx — Ci I e {x,j^ + Xj)dx. 
Then, C^ = 0. Thus, 

when |a| = 1. This also completes the proof of Theorem 15. II □ 



Corollary 5.1. Let d > 2. Suppose ()5.ip and (|5.3|) /lo/d. Further assume that there is a 
constant C such that 

|VVp <C(|Vt/p + l) (5.13) 

for all X e . Then, there is a constant A > such that 

d „ d „ 

V / \X,bj+Xjb^fdx>xY] / |X.,6jfdx (5.14) 
Jm^ ,j.^i Jm^ 

holds for any b = (61, 62, ■ ■ ■ ,bd) : M."^ M."^ E (^^(M'i))^ satisfying 



e 2 b{x)dx — 
e 2— (a;i6j — Xjbi)dx — 0, ^ < i ^ j < d. 
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Proof. Recall from (|5.6p and (|5.7p that one has the inequality 

d ^ 

V / \X,b,+Xjb,fdx 
> 2 / \Vbf + _ Iaf - ^ mv\)dx. (5.15) 



ij = l 

From the assumption (|5.13p . 



for any e > 0. This implies 



\WV\'>X{^\WV\'-Iav)-C, (5.16) 



for some constant < A < 4. Similarly, it holds that 



,1 



> (--Ce)|VFp-C, (5.17) 



for any e > 0. Then, combining (|5.16p and (|5.17p . it follows 

^\\/V\' -Uv-J2 mv\ > X{\\\7V\' - Uv) C, 

ij = l 

for some constant < A < 1. Plugging the above inequality into (|5.15p gives 

d 



V / \X,bj+X,b,\^dx 



ij = l 

>\[ \Vbf + \bfil\VVf -lAV)dx-Cx [ \bfdx. 
From the Korn-type inequality (|5.5p . it further follows 

V/ \X,b,+X,h\^dx>\ I |V6p + |6|2(l|VFp-iAF)da: 



aV / \X,bj\^dx. 



This proves (|5.14p and thus completes the proof of Corollarv l5.ll □ 

We conclude this paper with a remark that Proposition 15.11 Theorem 15.11 and Corollary 
15.11 can be described in terms of the probability measure as in [27], where 1^ 

satisfies the corresponding conditions. In fact, equivalently, (j5.4p become 



\ \^afd[i >x[ \a\^dn 
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for a with ad^ = 0. Similarly, (|5.5p and (|5.14[) are respectively equivalent with 

d ^ d ^ 



and 



^3 

for 6 = (61, 62, • ■ ■ I ^d) with 



d „ d „ 



bdfi = and / x x b dfi = 0. 
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